Nematic and chiral orders for planar spins on triangular lattice 
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We propose a variant of the antiferromagnetic XY model on the triangular lattice to study the 
interplay between the chiral and nematic orders in addition to the magnetic order. The model has 
a significant bi-quadratic interaction of the planar spins. When the bi-quadratic exchange energy 
dominates, a large temperature window is shown to exist over which the nematic and the chiral 
orders co-exist without the magnetic order, thus defining a chiral-nematic state. The phase diagram 
of the model and some of its critical properties are derived by means of the Monte Carlo simulation. 
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Nontrivial orders in frustrated magnets [1] are among 
the central issues in the field of condensed-matter physics. 
Besides the conventional magnetic order parameter of 
spin Si at a site z, there could appear various nontriv- 
ial orders such as vector [2, 3] and scalar [4, 5] chiral 
orders [6], and nematic order [7], which might lead to 
additional phase transitions distinct from the one driven 
by magnetic order. Even the ground state itself may be 
characterized solely by these nontrivial orders. This is- 
sue is now attracting revived interest from the viewpoint 
of nontrivial glass transition of spins [8] and multiferroic 
behaviors [9, 10], where the ferroelectricity is induced 
by the formation of vector spin chirality [11]. One im- 
portant aspect of this problem is the interplay between 
the various orders. Usually the nontrivial orders become 
long ranged when the magnetic order sets in. For exam- 
ple, the spiral spin order naturally implies the vector spin 
chiral order through (S^ x Sj) ~ (S^) x (Sj) on the neigh- 
boring sites. Therefore, the interesting issue is whether 
the nontrivial order can become long ranged in the ab- 
sence of the magnetic order. This issue has been studied 
theoretically [9], and experimentally in the quasi-one di- 
mensional frustrated magnet [12] where the chiral order 
appears above the magnetic phase transition. The next 
important question, we argue, is the interplay between 
the two nontrivial orders, e.g., chiral and nematic orders, 
which has not been fully addressed so far. 

To address this issue, we study a generalized classical 
XY spin model on a triangular lattice, 

H = JiY^ cos(%) + J2 XI cos(2%), (1) 

where Oij is the angle difference 9i — Oj between the near- 
est neighbors (ij). This model contains the usual frus- 
tration in the exchange interaction due to the triangular 
lattice geometry, together with the possible nematic or- 
der induced by the J2 term. The J2 = limit has been 



extensively studied, and it is believed to have two phase 
transitions at closely spaced critical temperatures [2, IS- 
IS]. The Kosterlitz-Thouless (KT) transition tempera- 
ture Tkt signaling the loss of (algebraic) magnetic or- 
der and the melting temperature of the staggered chiral- 
ity, T^, are extremely close, (T^ - Tkt)/T^ < 0.02 at 
J2 = 0, hampering the interpretation of the intermedi- 
ate, Tkt < T < phase as the chiral phase in which 
the chirality is ordered but the magnetism remains dis- 
ordered. Extension of the XY model to include large J2 
interaction was considered earlier in Refs. [16, 17], where 
the authors examined the phase diagram of Eq. (1) on 
the square lattice, which lacks frustration. In contrast, 
our model on the triangular lattice serves as a minimal 
model to study the two nontrivial orders, i.e., the chi- 
ral order induced by the geometric frustration, and the 
nematic order induced by the bi-quadratic interaction. 

A unique feature of the large J2 / Ji region of the model 
as noted in Refs. [16, 17] is the existence of an Ising 
phase transition associated with the vanishing string ten- 
sion between half-integer vortices in addition to the KT 
transition. This Ising phase transition turns out to cor- 
respond to the onset of the (algebraic) magnetic order. 
Being driven by Ji , the Ising transition temperature oc- 
curs at a much lower temperature than either the chiral 
or the nematic transition, which are both driven by J2. 
The result is the existence of a magnetism- free, chiral- 
nematic phase in the large J2/J1 part of our model. 

Phase diagram: The x — T phase diagram for Eq. 
(1) is shown in Fig. 1, where T is the temperature and 
X parameterizes the interaction as Ji = 1 — x, J2 = x. 
Detailed Monte Carlo (MC) calculations were performed 
with 5 X 10^ MC steps per run, on L x L lattice with 
L ranging from 15 to 60. Occasional checks were made 
on a larger lattice of up to L = 100 to ensure that no 
discernible changes in either the critical temperatures or 
the critical exponents are obtained from the larger size. 
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FIG. 1: (Color online) Phase diagram of the Ji — J2 model 
in Eq. (1) with Ji = 1 — x and J2 = x. Two closely spaced 
transition temperatures labeled by Tkt and separate the 
paramagnetic (PM) phase from the algebraically correlated 
phase at a lower temperature. aM, aN, and C stand for phases 
with algebraic correlations in (antiferro) magnetic and (anti- 
ferro)nematic order parameters, and the long-range correla- 
tions in the chirality order. A further transition from aN to 
aM occurs as an Ising transition for x > Xc with Xc ~ 0.7. 
All the symbols have a thickness in the temperature direction 
consistent with their statistical errors, (inset) The onset of 
chirality order at (red) takes place at temperatures close 
to, but slightly higher than the corresponding KT transition 
temperature Tkt (black) for all x. The Ising transition tem- 
perature Ti is not shown here for clarity. 

Typically, 10^ steps were discarded to reach equilibrium. 
An integer vortex-mediated KT transition marking 
the PM-aM boundary bifurcates into a half-integer 
vortex-mediated KT transition, marking the PM-aN 
boundary, plus an Ising transition [16] when x exceeds 
Xc ^ 0.7. The Ising transition in turn separates the aM 
from aN. For the whole range of the chiral transition 
temperature stays slightly above Tkt, with the 
possible exception at x = Xc where they may coincide. 

KT transition at Tkt- The determination of Tkt 
is made with the phase stiffness, also called the helicity 
modulus, appropriate for the J1 — J2 model 

-T^iiJi'^XijsmOij -^2J2'^Xijsm20ijf). (2) 

Here Xij = xi — Xj is the separation of the x- 
coordinate. The crossing of psiT) with the straight line 
(2/^) ( V^/2) ( Ji + 4 J2)T = (2/^) ( V3/2) (1 + 3x)T yields, 
for a given lattice size T, an estimate of the critical 
temperature Tkt (^) [14]. Extrapolation to T ^ 00 using 
polynomial fits as shown in the insets of Fig. 2 yields 
the estimate of Tkt- A more sophisticated method 
taking into account the logarithmic correction [18] yields 
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FIG. 2: (Color online) Phase stiffness ps{T) according to Eq. 
(2) for L = 15 - 60 and X = 0.5 and 0.9. The straight line 
is (2/7r) (a/3/2) (Ji +4J2)T. The crossing temperature of this 
line and Ps(T) for each L is shown in the inset along with the 
extrapolation to — 0. 

a similar answer [15]. 

Chirality transition at T^: It is customary to define 
the chirality x the directed sum of the bond current 
(sin6>ij)[13] following the relation (sin6>ij) ~ —dF/dAij. 
The free energy F is evaluated with respect to the mod- 
ified interaction cos Oij cos{Oij -\- Aij). A similar mod- 
ification of Eq. (1) results in the bond current 

J^j - Ji(sin(%)) + 2J2(sin(2%)). (3) 

This new definition is particularly effective as x 1, 
where the conventional definition ~ {sinOij) vanishes 
identically due to the Z2 symmetry. For each T^ was 
obtained from Binder cumulant analysis for the new def- 
inition of chirality based on Eq. (3). The conventional 
definition ( J2 = 0) gave an estimate of T^ which differs 
only in the third significant digit. Although our analysis 
showed T^ > Tkt for all x, we do not at present rule 
out the scenario in which T^ and Tkt merge at x = Xc, 
resulting in a multi-critical point there. If that happens, 
the second-order chirality transition may become weakly 
first-order. 

Earlier analysis [14] at x = identified the transition 
of X with the non-Ising critical exponents 1/u = 1.2, 
and /3/u = 0.12, j/u = 1.75. Figure 3 shows x 
and its variant, ip = {{x^) — {x)'^)/T^ in scaling 
form X = L-^/^f{tL^/^), i/j = L^^^ g {tL^^^) , with 
t = \T - T^\/T^, at X = 0.3 and x 0.8. Same 
exponents as for the x = case works well in scaling 
throughout the whole phase diagram. Appearance of 
the non-Ising exponents for J2 = have been explained 
in terms of an enhanced finite-size scaling effect at small 
sizes due to the screening length associated with the 
KT transition, in the cases of the square lattice [13] 
and triangular lattice [14]. Here it is equally possible 
that the true universality classes at T^ that of Ising 
transition. At any rate, the identification of the chirality 
transition T^ well above the magnetic transition for 
large J2 / Ji ratio is unequivocal and proves the existence 
of the magnetism-free, chiral-nematic phase in our model. 
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FIG. 3: (Color online) A scaling plot of chirality based on Eq. 
3) and its susceptibility for x = 0.3 and x = 0.8, for lattice 
sizes L = 15 — 60. The exponents used are those of x = 0[14]. 
The last row shows the behavior of the Binder cumulants at 
X = 0.3 and x = 0.8, respectively. 
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FIG. 4: (Color online) The size dependence of (a) the mag- 
netic (Ai) and (b) the nematic (A/") order parameters at 
X = 0.9 are shown on the log- log plot, (insets) The crit- 
ical exponent r]M{T) for M - l/L^^^^^ and m{T) for 



Magnetic and nematic orders: The low- 
temperature phase immediately below Tkt is either aM 
or aN, depending on whether x < Xc or x > Xc. The mag- 
netic and nematic correlations are examined on the basis 
of the order parameters, A4 = (3/i^^)| J^ie^ I' 
A/" = (3/L^)| Xl^^^ e^*^^ I, respectively, where the sum 
i e A spans the A sublattice sites. For Ti < T < Tkt, 
the magnetic order parameter is expected to lose its 
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FIG. 5: (a) A cartoon depicting the loss of head-tail order 
in going from aM to aN phase, (b) A snapshop of the chiral- 
nematic state at T = 0.2 for x = 0.9 where the Ising transition 
occurs at Tj = 0.177. Within the same sublattice the "body" 
of the arrows, not their tips, are seen to point in the same 
general direction. 



algebraic character and become short-ranged. Indeed, 
the size dependence of M as revealed hy M ^ 
for X = 0.9 has the exponents r]M{T) changing abruptly 
from ^ 1 above Tj to a small value below it (Fig. 4 

(a) ). The critical nature of the nematic order parameter 
AT at X > is seen in the continuous dependence of 
the exponent r]x^ ^ 1/L'^^^^\ as shown in Fig. 4 

(b) for X = 0.9. The T-dependent exponent r]^{T) 
continuously decreases as the temperature is lowered, 
even in the low-T magnetic phase T < Ti^ indicating 
that the nematic order remains critical in the whole 
temperature range < T < Tkt- A careful comparison 
of r]M{T) and r]^{T) for T below Tj revealed a relation 
rjjsfiT) ^ 4:7] m{T)^ in accord with the expectations of 
the spin wave analysis. 

Ising transition at Tj: A cartoon picture of the Ising 
transition is given in Fig. 5 (a), where it is described as 
the loss of local "head-tail" order. The choice of the order 
parameter for the transition is not unique and, to the 
best of our knowledge, has never been given an explicit 
expression. Here we choose to analyze the temperature 
dependence of 



J=(3/L2)^sgn(cos[^,-^,o]), 



(4) 



where Oio is the spin angle at some reference site iO of the 
A sublattice. As an Ising-like variable, sgn(cos[^i — Oio]) 
carries two allowed values ±1. In the aN phase, Oi and 
Oi -\- TT occur with equal probabilities, thus 1 = 0. An 
excellent data collapse in finite-size scaling was obtained 
with the 2D Ising critical exponents, /3 = 1/8, 7 = 1.75, 
and = 1 for both x = 0.8 and x = 0.9. To be exact, 
the orientation of Oi with regard to a reference angle 
Oio will be arbitrary as the separation i — iO tends to 
infinity in a truly thermodynamic system. Given the 
small exponent r]j\f{T) < 0.03 near T = Ti consistent 
with an extremely slow decay, however, one can argue 
that the only effective low-energy fluctuation is the 
TT-flip of the spin (which reverses the sign of cos[^i — ^^o]) 
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rather than the smah-angle fluctuations (which does not 
reverse the sign) for the practical system sizes considered 
in the MC simulation. As far as this is the case, our 
deflnition serves as a good measure of the Ising transition. 

A A A A 
A A A A 

FIG. 6: Eight possible magnetic patterns within the nemat- 
ically ordered phase, which also includes configurations with 
the global rotation of all the spins shown here. The corre- 
sponding chirality of each spin configuration is shown inside 
the triangle. 

Chiral-nematic phase: The central flnding of this 
work is the identiflcation of the chiral-nematic phase in 
the absence of any magnetic order. Although nematic 
order is algebraically ordered, the chirality, due to its 
discrete nature, can undergo a true long range ordering. 

The case for chiral-nematic order at large J2/J1 ratio 
can be made clearly if we consider the Ji — J2 model 
with the discrete angles Oi = 27rni/p, p = 6, and 
rii an integer between 1 and 6. The bi-quadratic J2- 
interaction turns into a 3-state planar model which is 
known to have a second-order transition (not KT transi- 
tion) into an ordered phase[19]. In our language, this is 
the paramagnetic-to-nematic transition. As the small Ji 
interaction is introduced, the six-fold spin model within 
the nematically ordered phase is governed by the effective 
interaction 

-(Ji/2)^(7,(7,, ai = ±l. (5) 

where the Ising variable denotes the two opposite ori- 
entations of the spin. Due to this residual interaction 
there will be an Ising phase transition at a temperature 
Ti ^ 3.641 X (Ji/2) ^ 1.82(1 - x) according to known 
results of the Ising model in two-dimensional triangular 
lattice. The linear decrease of T/ with x expected from 
the effective interaction as well as the absolute values of 
the critical temperatures are consistent with the phase 
diagram. Fig. 1. Within the nematic-ordered phase, 
there are eight spin conflgurations allowed for a triangle 
as shown in Fig. 6. The chirality for each conflguration 
reads 

X^k = i^i^j + ^j^k + Cr/eCri)/3, (6) 

using the Ising variables. For the downward triangle, the 
chirality is the opposite: xfjk ~ —{o'iCrj-\-ajak-\-(Jk(Ji)/3. 
Then, the net staggered chirality is given by 



X ^(cTiCr^- + djak + (TkCTi) ^aiaj. (7) 

(ij) 

The flnal expression, being proportional to the energy, is 
positive at any temperature T for a ferromagnetic Ising 
model given in Eq. (5). Therefore the chirality remains 
non-zero at temperature above Tj where magnetic order 
is lost, but the nematic order is long-ranged. It is possible 
that the chirality, due to its discrete nature, can survive 
the continuum limit p ^ oo and remain long-ranged even 
as the underlying nematic correlation is algebraically de- 
caying. 
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